Rotating Superconductors and the London Moment: 
Thermodynamics versus Microscopies 
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Comparing various microscopic theories of rotating superconductors to the conclusions of ther- 
modynamic considerations, we traced their marked difference to the question of how some thermo- 
dynamic quantities (the electrostatic and chemical potentials) are related to more microscopic ones: 
The electron's the work function, mean-field potential and Fermi energy - certainly a question of 
general import. 

After the correct identification is established, the relativistic correction for the London Moment 
is shown to vanish, with the obvious contribution from the Fermi velocity being compensated by 
other contributions such as electrostatics and interactions. 



I. INTRODUCTION 

The denning property of superconductors, the well 
known Meifiner Effect, is to expel an (undercritical) mag- 
netic field from its bulk. This property reverses itself 
when the superconductor is rotated, and a spontaneous 
magnetic field appears - again in the bulk, one or two 
penetration lengths away from the surface. This is usu- 
ref erred to as the London Moment!! Its magnitude 



B 



-[2m e (l + C)c/e]n, 



(1) 



where f2 is the rotational velocity, c the vacuum light 
velocity, m e the bare mass, and e < the charge, of 
the electron; finally, £ is a relativistically small correc- 
tion. All theories and experiments agree that £ is small, 
but strong disagreement exists with respect to its actual 
value and sign. The microscopiG-,tiieories, by AndersonB, 
BradyQ, Cabrera and coworkersBB, and BaymO take the 
main correction to be positive and given by the Fermi 
velocity, 



C w (v F /cf w 2 x 10" 



(2) 



The thermodynamic theory, on the other hand, finds (, 
to be negative, and very much smalleru: 



C = /i/c 2 « -10" 



(3) 



where jl is the chemical potential of the metal, the energy 
needed to add a unit mass to the superconducting solid. 
(Since the solid holds its atoms together, the chemical 
potential fi is necessarily a negative quantity.) The value 
of 10~ 10 makes the correction negligible for all conceiv- 
able purposes and renders Eq([l]), taking the values of e 
and c as given, a very precise expression for m e - about 
three orders of magnitude more precise than any present 
experimental technique to determine the electron mass 
directlyta. 



Generally speaking, although the thermodynamic, 
macroscopic theory lacks details and is incapable of an- 
swering many quantitative questions, it is nevertheless 
a rigorous theory, and uniquely appropriate for under- 
standing the London moment. This is because (i) the 
London Moment is an equilibrium phenomenon, (ii) the 
measurement concerns a quotient between two macro- 
scopic fields, B and S7, and most importantly, (iii) the 
reason for the microscopic parameters of m c and e to ap- 
pear on the macroscopic level is due to simple symmetry 
principles, rather than some detailed properties of the 
interaction. 

Comparing both types of theories in detail, we find 
their difference to be easy to bridge, if the relations be- 
tween some thermodynamic and microscopic quantities 
had been clear beforehand. More specifically, it is the in- 
correct identification of electrostatic and chemical poten- 
tials on one hand, the work function, mean-field potential 
and Fermi energy of the electron on the other that has 
led to the above discrepancy. The correct identification 
of these quantities is of rather general interest and tran- 
scends the understanding of the London Moment alone. 
The proper identification is the most important result of 
the present paper. 

Another problem clouding the understanding of the 
London Moment is the disagreement between experiment 
and all theories. Although this remains a point we do not 
understand, we took a step towards its clarification by 
achieving understanding of the following point: The best 
present experiment measures the flux_^ A • ds = J B • da, 
outside the rotating superconductor!!!]. Two effects con- 
tribute the this flux, the London Moment, and the so- 
called double layer. The latter is a result of the fact that 
a metal may be conceived of as an electrostatic poten- 
tial of the square-well form, with the depth A$. As 
the metal rotates, this discontinuity in the rest-frame 
electrostatic potential produces a discontinuity in the 
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FIG. 1: Direct measurement of the London field Bl- 



laboratory-frame vector potential, 



AA = -(SI x r)A$/c. 



(4) 



And this also contributes to the flux § A • ds, with a 
contribution much larger than the relativistic correction 
of the London Moment. It has been generally suspected 
that both contributions are experimentally inseparable, 
that they are always measured in conjunction. (A double 
layer of opposite charge - or one layer of dipoles - pro- 
duces a discontinuity in the electrostatic potential, hence 
the name.) 

We disagree. In this paper, we shall discuss both ef- 
fects separately, as only the first is universal - the second 
varies with materials, and is quite independent from su- 
perconductivity. Besides, we believe that these effects 
may in principle be measured independently. Taking the 
superconductor to be a cylinder rotating around its axis, 
the London Moment may be measured immediately out- 
side its top or bottom surface, where any magnetic field 
B is continuous, see Fig 1. (It is not continuous across 
the cylinder surface, where the field drops to zero within 
a few penetration lengths, due to the presence of per- 
sistent currents.) The flux from the double layer alone 
may be seen above T^jOr in the experiment proposed in 
chapter 0. In chapter |l[ we shall present the simple yet 
stringent macroscopic calculation leading to the London 
Moment. This is followed in chapter jTj by a compar- 
ison to the microscopic results, which are brought into 
agreement with the macroscopic ones by the appropri- 
ate identification of the above mentioned thermodynamic 
and microscopic quantities. In chapter IV, we consider 
the effect of the double layer, the contribution of which 
may be measured in the experiment described in chap- 
ter jy|. 



II. THE MACROSCOPIC APPROACH 
A. Electrostatic and Chemical Potential 

Given any neutral, macroscopic system, we may either 
change its numbers of electrons and ions, 2V_ and JV+, or 
its mass and charge, M and Q. They are related as 



M = miiV + + to c AL 
Q = \e\(N + -N_). 



(5) 
(6) 



(mi and m e denote the bare mass of the ion and electron, 
respectively, and e < is the elementary charge.) Fre- 
quently, the number of atoms N a = Mj (mi + m e ) is used 
instead of the mass M. The energy of these changes are 
given by the respective chemical potentials, 



dE 



where employing Eqs 
to be related as 



H+dN+ + Li-dN- (7) 
fidM + A$dQ (8) 
MadA^a + A$dQ, (9) 

) , we see the chemical potentials 



(10) 

(11) 
(12) 



jU+ = TOj/i — eA$, 
/i_ = m /i + eA$, 
/i = /i a /(m; + m ). 



Because //_, //+, /i a , [i and A<I> respectively denote the 
energy needed to bring an electron, an ion, an atom, unit 
mass and unit charge from infinity to inside the system, 
these are fixed quantities, and we are not at liberty to 
alter them by an additive constant. We especially have 



A$ = $(r) - $(cx)), 



(13) 



where $ is the electrostatic potential, which is fixed 
only up to an arbitrary gauge transformation, and which 
forms a 4- vector with the vector potential A. Note that 
the energy and chemical potentials employed in this pa- 
per include the rest mass, because otherwise the relativis- 
tic transformation formulas of the following sections will 
not work. Hence we define the non-relativistic chemical 
potentials with tilde, 



= m c c 2 + ji-, 
fi = c 2 + /i, 
[i a = (mi + m c )c 2 + fi a 



(14) 
(15) 
(16) 
(17) 



Usually, both fi a and /2_ are of the order of a few eV, 
this being the scale of atomic physics, cf an early model- 
calculationO. Since m c /(rrii + m c ) ~ 10~ 5 , Eq ( |TT| ) and 
( |l2l ) imply, to great accuracy, 



/i_ = eA$, 



it.h 



fj.-/m c 



10 



fi : 



10- 10 c 2 . 



(18) 



(19) 
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B. Josephson Equation and Superfluid Velocity 



The Josephson equation is usually given as 
(&/2)£ + /i- = 0, 



(20) 



though this equation is not gauge invariant. In view of 
Eq ( |ll| ) and (fL3l), the correct form for the Josephson 
equation is clearly 



(h/2)(p + m e /i + e$ = 0, 



(21) 



implying that Eq ( |20| ) is valid only if $(oo) = 0, a special 
gau ge c hoice. (Similarly, taking (p = when considering 
Eq ( pi| ) also constitutes a gauge choice.) Note that the 
form of Eq (|2l]) does presume rest frame of the crystal, 
and equilibrium. The superfluid velocity is also defined 
as a gauge invariant quantity 



1 

m e 2 



:a). 



Its equation of motion reads 

v s + V/i = eE/m e , 



(22) 



(23) 



obtained by inserting Eq (|22j) in (21) and employing E = 
— (V$ + A/c). Although this pleasingly simple equation, 
independent from any coefficients, may be used to justify 
the definitions of Eqs (21 2^), we must be aware that 
v s as defined need not transform as a velocity under a 
Galilcan-Lorentz boost. The transformation behavior of 
v s is given by the Josephson equation (|2l|): Define a 4- 
velocity 



Ti dtp 



2m e dx a 

h 



m c c 
((p/c, V(p) 



-A a 



2m e " ' ' m c c 

and use Eq ( ^i"|) to yield 

u a = (-///c, v s ). 
Now, because u a transforms as 

u' a = (uq — u-v/c,u - u v/c) 



(24) 
(25) 

(26) 



under a boost of d, to linear order in v/c, so does the 
4- vector (—^i/c,v s ), leading to 



,s\l 



/' - 

V s 



{n/c 2 )v, 



implying that it is the quantity 



( M /c 2 )-V = ^/(i + c) 



(28) 
(29) 



(30) 



that transforms as a velocity, (v ts )' = v ts + v. We shall 



refer to v ts as the true superfluid velocity 
with 



C = A/c 2 = (fi- - eA$)/m c 2 



10 



-10 



although 
, (31) 



cf Eqs (0jT6|JT|), the difference is very small indeed. 
Defining an effective mass to* by 



TO* = TO C (1 + 



1 6 A N 

m* 2 c 



(32) 
(33) 



we shall find in the next section that this is the mass 
appearing in the formula for the London Moment, Eq (Q). 



C. Thermo- and Hydrodynamics 

In a general inertial frame, the thermodynamics of a 
superconducting solid is given as 

de = T ds + fidg + v n ■ dg + atjdVjUi 

+E • dD + H • dB + j s • dv s . (34) 

Notation and explanation: s is the entropy density, and 
T the temperature, /i dg and E° • dD are the respective 
local expressions for /idM and A$dQ of Eq (||), where 
E ■ dD involves a partial integration, and is manifestly 
gauge invariant; H ■ dB is the magnetic counter term. 
j s • dv s characterizes the broken phase invariance of su- 
perconductivity, as (TijdS/jUi the broken translational in- 
variance of solids - with Ui the displacement vector. In 
a general inertial frame, the total, conserved momentum 
density gr tot is also a thermodynamic, variable, though we 
employ g = g tot — D x B insteadt 2 !, with v n = de/dg 
being the equilibrium velocity of crystal points, atoms 
and normal electrons. In the present case of interest, 



O x r. 



(35) 



Being a conjugate variable, j s = de/dv s is a function of 
those two thermodynamic variables also odd under time 
reversal, v s and v n . In an expansion, to linear order of 
the variables, we have 



(36) 



where and otjk are two expansion coefficients, while 
the prefactor (c 2 /fi) simply renormalizes p\ -. A Maxwell 



(27) relation then implies 




(37) 



Confining ourselves to the local rest frame and disregard- 
ing dissipative terms (then s, = 0), the hydrodynamic 
set of equations is given by the Josephson equation (|2^), 
the Maxwell equations, and the conservation laws for en- 
ergy and mass 

i = -V Q, Q + V-i p = 0, (38) 
B = -cVxE, D = cVxH-j e , (39) 

where Q, j p and j c are as yet unknown. Inserting these 
expressions into the temporal derivative of Eq (|3~1), e = 



4 



fig +ED +H-B+ j s -V s , and insisting that all equations 
are satisfied simultaneously, we find 



j P = js, jo = ej s /m c 
Q = jsM + cE x H. 



(40) 
(41) 



Clearly, j s = de/dv s has the significance of being the 
persistent mass and electric current in the rest frame. 
The relation 



i P = ej s /m e , 



(42) 



is a necessarily one, as the transfer of one electron is 
coupled to the transfer of the bare values of e and m e . 
Starting from this relation and tracing it back, we would 
have found that the charge e in Eq ( p2| ) must indeed 
be the bare charge, and that v s must indeed be gauge 
invariant. Because of the symmetry of the energy stress 
4-tensor, g tot — Q/c 2 , we have 



isH/c 2 



(43) 



for E, D 
implies 



(and still in the local rest frame). This 




Ex 



or in combination with Eq (|3 

it =pi i [(c 2 /n)v°-v»}=pt j [v 



(44) 



ts 
1 3 L "J 



Inserting this expression into Eq ( |38| ) with D = 0, taking 
another curl on both sides, and denoting the matrix p|- 
as /5 s , we finally obtain 



V x 



— (p")- 1 V x h| = V x {v ts - v n ) 
e 

m e (l + Qc 

This is the equation that accounts for the equilib- 
rium behavior of the magnetic field in superconductors: 
(m e c/e)(p s ) -1 yields the square of the inverse penetra- 
tion depth, a quantity that depends on the crystal di- 
rection, while the vanishing of the right hand expression 
gives the bulk value of the magnetic field, Eq (Q). Be- 



cause of Eq (31), we may generally neglect the factor 

i + C 



III. THE SQUARE- WELL POTENTIAL 

Consider a square- well potential, Fig 2, a popular- 
model for metals. Taking the outside value of the en- 
ergy as zero, the depth of the well is eV(< 0), with V 
the mean-field potential for electrons; filling the potential 
up to the Fermi energy e F (> 0), the gap still remaining 



may be identified as the work function of the electron 
W(< 0). Together, they satisfy 



W -e F = eV. 



(47) 



When making contact with thermodynamics, there are 
no doubts that the identification 



W =/L 



(48) 



holds, because the physical significance of W and the 
energy needed to take out an electron, and the energy 
gained when putting one in, respectively, is the same. 
The other two scales, the Fermi energy e F and the mean- 
field potential V, are model-dependent quantities, with 
a large measure of arbitrariness attached to them. So 
we should not expect them to be directly measurable, 
or connected to thermodynamic quantities in a simple 
fashion. Nevertheless, all mentioned authors employing 
the rajiSiMscppic approach to understand the London Mo- 
menti2o&Q0 adopt the identification 



V = A$, 



(49) 



with some apparently plausible consequences. First, with 
Eqs <M) and @, we have m* = m e (l + C) with 



C(m c c 2 ) = Jx_ - eA$ = W - eV = e F . 



(50) 



This makes the mass correction positive, and rather tan- 
gible: With the estimates W « -4eV, eV « -96eV, 
and e F ~ 92 eV (by averaging the electron wave func- 
.*], (45) tions over_the Fermi surface), the value ( a 1.8 x 10~ 4 
was founcH for the mass correction - as compared to £ ~ 
— 10~ 10 of the last section. Moreover, with e F — ^m e Vp, 
Eq ( [so]) delivers the simple, kinematic-relativistic inter- 
pretation for the mass correction, 



C = \{v F /c) 2 w 1.8 x 10" 



(51) 



So why do we claim that Eqs ( |49] , |50] , |51| ) are in error? 
The identification of Eq ([i9|) is made by taking the 
macroscopic electrostatic potential A$ as the potential 
felt by superconducting electrons from "all charge dis- 
tributions in the metal, such as surface dipoles..., the 
screening hole..., and charge inhomogen&ities associated 
with atomic cores and valence electrons' u. (Sentences in 
italic are quotes, here and below.) Clearly, the one os- 
tentatiously lacking in this list is the contribution from 
the other band electrons (including the superconducting 
ones), to the potential, which brings the potential energy, 
gained by the last electron to be added to the system, to 



W = fi- = eA$, 



(52) 



cf Eq (|18|) . This is the thermodynamic definition of the 
potential A$, see the discussion of the last section, and 
we are not at liberty to alter it. In fact, as already men- 
tioned, it would have been highly surprising for the di- 
rectly measurable macroscopic electrostatic potential $ 
that forms a 4-vector with the vector potential A to be 
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Metal Surface 



e*3> t . . =0 

outside 



w 



5<D 



inside 



Mean-field 
Potential ■ > 



eV 



FIG. 2: The relations between the macroscopic electrostatic 
potential e<3>, the work function W , the Fermi energy ef, and 
the mean-field potential eV, close to a metal surface. 



simply related to the model-dependent mean-field poten- 
tial V. 

This still leaves the kinetic-relativistic interpretation 
of the Fermi velocity, Eq (|5l]), intact. As it is a rather 
popular preconception, and employed by two more-cecent 
papers without the detour via the Fermi energytLij, we 
need to emphasize here that it is really only correct for 
free particles, and does not work for strongly interacting 
systems. 

To understand this, consider an atom, consisting of an 
electron and an ion - the same way the metal consists of 
band electrons and the rest. We may write the energy of 
the atom as E = (m* + m* c )c 2 , or 



(53) 



where the superscript * denotes the effective mass includ- 
ing the mass defect, while W the binding energy between 
the electron and the ion, the analogue of the work func- 
tion above. Note that the binding energy W includes all 
contributions, especially the kinetic energy of the elec- 
tron. Microscopically, we may divide W between to* and 
to* arbitrarily. But if the division is in proportion to the 
rest mass, a natural (if not obvious) procedure, the mass 
defect of the electron is 



W 



(54) 



in agreement with the result of the last section, see 
Eqs PJl7p|). 



Baym considered the problem of mass correction 
within the Landau theory of Fermi liquid in his workQ. 
The theory is beyond the usual mean-field approach, and 
accounts for many-body effects. Because his arguments 
and calculations are complicated, we shall only quote his 
final expression, 



(to* - m c )c 2 = e kin + e lnt - eV c + 



(55) 



where £/-j n is the kinetic energy of a Fermi electron, sim- 
ilar to Ep above, while the last two terms are new: V c is 
given by "the average of the electrostatic potential within 
a given unit cell (eg the Wigner-Seitz cell), arising from 
the charges within the cell. " And ei n t is "the interaction 
energy of an electron at the Fermi surface, containing 
electrostatic as well as exchange contributions (but ex- 
cluding the boundary dipole layer)." The value of eV c 
was given as — 19.8ey, while the value of Si n t remains 
unknown. Because e^t not only includes many-body in- 
teractions, but also electrostatic contribution, it may well 
be an appreciable term, canceling other contributions to 
yield virtually vanishing total correction, in agreement 
with the thermodynamic result. Baym also points out 
that corrections from the presence of the lattice are to 
be expected, because the Landau Fermi liquid theory is 
confined to translationally invariant systems, of which 
metal electrons are not one. Only a generalized Landau 
theory that also considers the ions would indeed be deal- 
ing with a system that is, in its totality, translationally 
invariant. This should result in more generally valid re- 
lations, Replacing relations such as g = px_ /(m e c as 
given inO, with g = \p a /(m e + mi)c 2 ]j s: or Eq (ph. 



IV. THE BOUNDARY DOUBLE LAYER 

Because of the presence of surface dipoles, the elec- 
trostatic potential is discontinuous at the surface of the 
superconductor, A4> ^ 0. If the metal is in rotation, the 
vector potential will also be discontinuous, contributing 
to the total flux of the system in the laboratory frame. 
Since the magnitude of this is ~ A<£>, taking eA<£> either 
as ^ = W w -4eV, Eq (HHf), or as eV ~ -96eV, 
Eq (^9|), clearly makes a big difference. As argued at 
length in the last two sections, we believe the first is cor- 
rect, as does the bookEJ by Landau and Lifshitz, who 
unequivocally pronounced the equality between the work 
function and the potential's discontinuity, W — eA$, see 
§23. 

The jump of the vector potential follows from the 
Lorentz transformation property of the electromagnetic 
4-potential A a . Consider a frame co- moving with the 
metal. If A' denotes the electric potential outside the 
metal, the potential inside the metal will be A' Q — A$. 
So the 4-potential A a in the local rest frame is (A Q , A ) 
outside the metal, and (A — A<£>, A') inside it. In the 
laboratory frame, the 4-potential A a becomes, in linear 



6 



order of v/c, 

(A - v • A /c,A - A v/c) 

outside the metal, and 

(A' - A$ -v • A /c, A - A' v/c+ A$d/c) 

inside the metal, where v is the velocity of the boundary. 
So the discontinuity in the vector potential is 

A A = A out ~ A mt = -(A$/c)w (56) 

For a metal cylinder uniformly rotating about the z-axis 
e z , we have v = Qe z x _R, with R the position vec- 
tor of the boundary. So the magnetic flux of a non- 
superconducting, rotating metal is 

<j)AA-ds= -2nR 2 A<f>n/c. (57) 

Since the experimentS measures both this boundary flux 
and the London field simultaneously, it is convenient to 
introduce an observed mass m bs (if the effect due,to the 
penetration depth is negligible, as is the case inEj). We 
write it as 

m ohs = (1 + C)(l - a)m c = (1 - a)m* c (58) 

where a accounts for the flux from the boundary. .double 
layer. In the cited microscopic considerations,&Q'ElJ it is 
the Fermi energy, ep w 92eV, that enters £, and the 
mean-field potential, eV « — 96eV, that enters a, though 
in opposite directions. Combined, the observed mass is 
corrected only by the work function W = eV+ £j? « 
— 4eV, cf Eq (f|7j), leading to a total correction of about 
8 x 1CP 6 . Thermodynamically, although £ m — 10~ 10 is 
negligible, but the work function W does enters a, so the 
measured effect is again given by a value around — 4eV. 

Since the same value for m b s is predicted by all theo- 
ries, one may conjecture that the values of ( or a individ- 
ually are unimportant, because any experiment can only 
observe {m t, s /m e ) — 1 = £ — a. This does not seem right 
to us. First, both the strength of surface double layer of 
a metal and the London field are well defined physical 
quantities. And the parameters £ and a are unambigu- 
ously related to the chemical potential and the work func- 
tion, respectively. Second, both effects are (at least in 
principle) independently measurable. The London field 
could be observed directly by measuring the flux near the 
top-center of a rotating superconductor (see Fig. 1), and 
the flux from the double layer only may be observed in 
the normal state, or even in the superconducting state, 
by the method given in the next chapter. 

V. A PROPOSED EXPERIMENT 

Consider two concentric, co-rotating, hollow cylinders, 
made of the same superconducting metal. The electric 




FIG. 3: Apparatus for measuring the boundary flux in the 
superconducting state. Superconducting metal is depicted by 
shaded area. 

current J flowing through the inner cylinder is measured 
by a SQUID, see Fig 3. The cylinders rotate with the 
angular velocity il. We calculate the magnetic field and 
the current in the apparatus. For simplicity we assume 
that the cylinders are of infinite height, and no surface 
charges (or electric field) are present. n 

This geometry was firstly analyzed by Bradyu, who 
estimated the effect of penetration depths A by sim- 
ply assuming constant magnetic fields in the regimes: 
< r < R\ + A and R 2 — A < r < R3 + A, in other words, 
by assuming that the actual internal spaces, < R\ and 
i?2 < r < i?3, are increases by A at the superconduct- 
ing boundaries. Moreover, because the magnitudes of 
the field was estimated by dividing the magnetic flux by 
the respective area, the contributions from the boundary 
layer have also been neglected. 

The current J in the inner cylinder as determined by 
Brady is proportional to the observed mass m b s and the 
difference of the angular velocities of the two cylinders. 
This was given as a method to measure the London field. 
In our case, both angular velocities are the same, so J 
would vanish in his approximation. Our more detailed 
calculation, however, shows that J is finite and observ- 
able if the distance between the two cylinders, R3 — R2, is 
small enough. In fact, it is proportional to aQ, rendering 
a, or the flux from the boundary layer, measurable. 

The magnetic fields and the current density inside 
the cylinders can be obtained by solving Eq ( |46| ) in 
the cylindrical symmetry. Considering the isotropic case 
Pij — Qs&ij with g s the superfluid density, we write the 
vector fields as: 

B = Be z , A = A(e z xr/r), 

f = f(e z x r/r). (59) 

The general solution for H = B is 

B = 2-/*n + dl + C 2 K . (60) 
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7 



-m*c/e; Ci_5 are integration con- and at R 3 : 



Notations: 

stants; A = (m e c/ |e|)\/ (1 + 0/ Qs is the penetration 
depth; In,Kn (N — 0,1) are Bessel functions with the 
arguments r /A: 0o = —hc/2e is the quantum fluxoid. 

Inserting ( |60| ) into Eqs B = V x A, j° = cV x B, and 
noting the quantization property § Vy? • ds = 27rn for the 
phase in Eq ((33|), we obtain the magnetic potential and 
current density: 



A = 7*fir -n<p /r + \{C\h 
j e = -(c/AXd/i-C^). 



(61) 
(62) 



The integer n denotes quantum state of the cylinder. The 
induction B is constant outside the rings, while the po- 
tential is of the form 



A = Br/2 + C/r 



(63) 



with C again a constant. Using the above solutions, it 
can be easily seen that the field of our apparatus is 



B = 



B 2 ; 
2j*n 



Cih 

C3I0 

and the potential is 



C 2 K ; 
C 4 K ; 



r < 

i?i < r < R 2 , 
R 2 <r < R 3 , 
R 3 <r < R 4 , 



(64) 



B ir /2; 
7 *Or - mfo/r + A(Ci/i 



A= < 



r < R t , 
-i-Caifi); 
i?i < r < R 2 , 
B 2 r/2 + C 5 /r; R 2 <r<R 3 , 
j*nr - n 2 <fi /r + \{C 3 h ~ C 4 K X )- 

R3 <r < R 4 . 



(65) 



Note that ri\, n 2 may be different, because the two cylin- 
ders can be in different quantum states. The current per 
unit length in the inner cylinder is 

J = / fdr 

J Ri 

= -c[Ci(/ 2 - hi) + C 2 (K 02 - K 01 )}. (66) 

Here and below we will use the notations 

Inm = In{Rm /A), Knm = Kn(Rm /A), 

where the integers N = Q,\, M = 1 — 4. The unknown 
constants Ci_5, B\,B 2 are determined from the bound- 
ary conditions. At according to Eqs (|64],|65|), we have 



B x = 2 7 *0 + d/01 + C 2 K 01 , 
BtRt/2 = (7* - A$/c) 
-n^o/Rx + A(CiJ n - C 2 K U ), 



(67) 



(68) 

where Eq ( |67| ) shows that the induction varies continu- 
ously at the surface. Eq ( |68|) is obtained from the jump 
of the potential Eq ( p6| ) . Similarly, we have at R 2 : 

B 2 = 2 7 *fi + C1/02 + C 2 K 02 , (69) 
B 2 i? 2 /2 + C 5 /R 2 = (7* - A$/c) flR 2 

-n 1( f> /R 2 + X(dl 12 - C 2 K 12 ), (70) 



B 2 = 2~f*n + C 3 / 03 + C 4 Ao 3 , (71) 
B 2 R 3 /2 + C b /R 3 = (7* - A$/c) Oi? 3 

-n 2 0o/i?3 + A(C 3 Ji 3 - C 4 A 13 ). (72) 

The induction is zero beyond R4. (We consider the case 
of no externally applied field, ie all the sources of the field 
are from supercurrents in the rings.) We have then 



2 7 *0 + C 3 I 04 + C 4 K 



0. 



(73) 



Solving Eq (B7-73) for Ci_5, B 4l B 2 and inserting them 



into Eq (|66|), we obtain the current. The result can be 
written as 

J = [/sni + U(n 2 - m)]J + c 7 *(/ 2 - fia)Q, (74) 



ith 



3A$ 



mzc* 



Jo 



Tic 2 
2eR\ ' 



(75) 



Here the factors /1-4 depend only on the penetration 
depth and geometry of the system. For observing the 
boundary effect, it is convenient to choose the geometry 
such that I/2 1 << \cxfi\. One example satisfying the 
requirement is: R\ — 2.5 cm, R 2 — Ri = 1 /im, R 3 — R 2 = 
40 nm, i?4 — i?3 = 1 /im, A = 40 nm. We have obtained, 
by the numerical computation, the values 

fx = 1.33333, f 2 = 1.852 x 10 -11 , 
f 3 = -1.9999968, / 4 = 2.083 x 10 7 . 

Because f 4 is very large, it is not probable for the two 
rings to lay in different quantum states. So within a good 
approximation, (and also neglecting £ <C 1) we have 



J=-2J ni-/i(A$)n, 



(76) 



with /1 = 4/3 for the geometry considered. This shows 
that the strength of the surface double layer A$ can be 
measured by measuring the change of the current with 
the rotation. Note that the factor f\ decreases with in- 
creasing distance R 3 — R 2 between the two cylinders, see 
Fig 4. When the distance is large, fi —> 0, the last term in 
(|74|) is negligible, and we return to the result by Bradyu. 



VI. CONCLUSIONS 

The conclusion of this paper is that the relativistic cor- 
rection to the London field is not the result of the Fermi 
velocity. Instead, it is given by the chemical potential fi 
of the metal, which quantifies the complicated interaction 
among all the particles, including that between electrons 
and ions. Because the interaction energy is typically sev- 
eral eV per atom, while the mass of an atom is around 
10 4 MeV, the relativistic correction is tiny, of order 10~ 10 , 



8 



fi 




FIG. 4: Variation of the factor /i with the separation distance 
of the two rings. 



and beyond the scope of any present experimental tech- 
niques. When previous, microscopic theories considered 
the London Moment, they considered the kinetic con- 
tribution to the correction £ of the electron's mass, but 
neglected other contributions, especially from the inter- 
action between the electrons and the lattice. In addition, 
the discontinuity of the macroscopic electrostatic poten- 
tial at the metal surface was incorrectly taken to be the 
mean field potential, or the sum of the kinetic energy and 
the work function. As soon as these errors are revised, 
£ is found to essentially vanish. Unfortunatelvp-the rea- 
son for the discrepancy between the experimented and all 
theories remains unclear. 
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